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The Thompson cross section for scattering of electromagnetic waves by a free electron in an 
expanding universe is derived here. The equations of motion of the electron are obtained from the 
Einstein-Maxwell field equations of general relativity using the Einstein-Infeld-Hoffmann surface 
integral method. These integrals are evaluated approximately by perturbing off an Einstein-deSitter 
cosmological field. It is found that the Thompson cross section varies as the inverse square of the 
cosmic scale factor R{t). 
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. . . The question of whether or not the constants of nature actually vary on a cosmological time scale was first inves- 
C ' tigated seriously by Dirac. j^] To satisfy his Large Number hypothesis, Dirac assumed that it was the gravitational 
. constant G that varied on a cosmological time scale. Later [|] he developed this idea further by introducing the 
' notion of gravitational and atomic time and assumed that these two times were not synchronous with each other. |^ 
However, in a recent work |^ (referred to hereafter as Al), I showed how one could construct model gravitational 
and electromagnetic clocks whose dynamics in an expanding universe follow directly from the Einstein-Maxwell field 
] equations of general relativity using the Einstein-Infeld-HofFmann (EIH) H] surface integral method without other 
assumptions. It was found that, for a class of clocks for which Tl ^ Tc ^ Th where is the light travel time across 
the clock, Tc is the clock period and Th is the Hubble time, that, to an accuracy of 0{e^6) where e = Tl/Tc and 
' S = Tl/Th, the two clocks remain synchronous with each other and measure the same cosmic time t when the cosmic 
, ■ gravitational field has the form 

g^,= dmg{l,~R\-R\-R^). (1) 

m ! . 

, In this paper I address a different problem associated with physical processes in an expanding universe, namely 
^\ • Thompson scattering. I will show that the Thompson scattering cross section is in fact not a constant but varies with 

[ time by a factor proportional to 1/R^{t). 
C ■ I At first glance it might appear to the reader paradoxical that gravitational effects associated with the expanding 
' universe could have any effect on what is seemingly a purely electrodynamic process. But in fact Thompson scattering 
O . is not a purely electromagnetic process. The inertial force term in the equations of motion for the scattering electron 
O^' is of gravitational origin as EIH first showed and, in the case of an expanding universe, involves the scale factor R(t). 
, As in Al the exact dependence of the inertial force on R{t) as well as the form of the electromagnetic force due to the 
incident wave can be determined from the EIH surface integrals by perturbing the gravitational and electromagnetic 
" fields off the Einstein-deSitter field as in Al. However, since one is dealing here with radiation fields it turns out to 
• ■ be more convenient to employ conformal coordinates (r, x, y, z) in which this field has the form 

g^, = R\T)Aiag{l, ^1,-1,-1) (2) 

and where r and t are related by dt = R{T)dT. 

The EIH surface integrals can most easily be derived from the Landau and Lifshitz Q form of the field equations 
for g^i, and the electromagnetic field F^^, and have the form (in what follows I will use units in which G — c — 1 , latin 
indices run from 1 to 3, greek indices run from to 3, and I employ both the Einstein summation convention and the 
comma notation to denote partial derivatives) 



where 



and 



iDTT 



and 



pf,p\ 



0. 



(5) 



In these equations g — det((ji^^), g^'^ 



and T'^'^ is the electromagnetic energy-stress tensor given by 



-gF^^ , l^j^ is the Landau-Lifshitz energy-stress pseudotensor 
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Note that, in our expression for O^"^ above, there is no matter contribution from the sources since they are assumed to 
be compact and to vanish on and outside the EIH surfaces. This feature of the EIH method thereby avoids having to 
make specific assumptions about the form of the matter energy-stress tensor or the need to introduce singular source 
terms. 

Because of the antisymmetry of F''" and U'^'^'' in their indices, integration of Eq. (||) over a closed 2-surface in a 
t = const, hypersurface gives 

^(C/^''°,o-e'"')n^dS' = (7) 

where rij. is a unit surface normal. In a similar way one gets from Eq. (j^) the result 

r\nnrdS = Q. (8) 

It is these last two equations that yield equations of motion for the sources of the gravitational and electromagnetic 
fields. 

When the surfaces over which the integrals in Eqs. (Q) and (^) are taken surround a source, the requirement that the 
surface independent contributions vanish yield these equations. (The surface dependent terms will in all cases vanish 
either identically or as a consequence of the field equations. Q) These equations will be used to derive equations of 
motion for a compact charge in the presence of an incident electromagnetic field whose characteristic length scale is 
large compared to the size of the EIH surfaces needed to enclose the charge. The near electromagnetic field produced 
by the motion of this charge is then matched to a far radiation field. The total flux of this scattered field is finally 
divided by the incident flux to yield an expression for the Thompson scattering cross section. 

To evaluate the fields appearing in Eqs. (zh and (^) I follow the methods used in Al except that here the gravitational 
field is perturbed off the background field rather than the field . For convenience the total field is written as 

and g^'^as well as the electromagnetic four-potential is expanded in a double series in e and S. These fields in 
addition are assumed to depend on the conformal time r through their dependence on er and St while i? is a function 
only of 6t. (In higher orders of approximation they will of course depend on higher order multiple times.) In addition, 
charges and masses are scaled so as to be O(e^). In what follows we will need an accuracy of and eS since we are 
only concerned here with the modifications in the Newtonian dynamics in the background cosmological field. 
The lowest order correction to the gravitational field is taken to be ||] 

5°" = 1 + £^h (10) 

and satisfies 

V'^h = 0. (11) 

For compact spherical sources (and Schwarzschild black holes) the solution in the weak-field zone on and outside the 
EIH surfaces has the form 

h^iy^ (12) 

where the index A labels the sources in the system and the sum is over all A. The m^iare as yet to be determined 
functions of er and St and r^i = x — xa where the are the coordinates of the Ath particle and are also functions 
of £T and St. When the surface integrals in Eq. ^ with /i = are evaluated using this field one finds that 

derffiA — and RdsTfriA — —dsTRffiA (13) 

so that, to this order of accuracy, 

ruA 

"^A = (14) 
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where the tua are constants. In a Hke manner one constructs the lowest order contribution to the scalar potential 

^° = (15) 

with (j) satisfying 

72 



V^0 = O. (16) 



When the spherically symmetric solution 



= E- (17) 



is substituted into the surface integral (9) one obtains the result that 

d^rQA = dsrQA = (18) 

so that 

QA = qA (19) 

where the are constants. 

To derive equations of motion from the surface integrals (^ we need the lowest order corrections to g"'' which are 
0(£^) and 0{e'^6) so we set 

=£^hl+e^6h^s. (20) 

There will of course be additional corrections in higher orders of approximation that are small compared to the 
first correction but large compared to the second, e.g. post-Newtonian corrections of order e^, but here we are only 
interested in the first order effects of the expanding universe on Newtonian physics. As in Al and are determined 
to be 

hi = 4:y2—XA,er and /i^ = 4V— a;^,5^ . (21) 
^ rA ^ TA 

when the gauge conditions 

hl,r +derh = and h's,r + ffiAdsr (J-^ = (22) 

are employed. 

Since my purpose here is to examine Thompson scattering in an expanding universe I will, in what follows, confine 
my attention to a single charged source and take the electromagnetic field used to evaluate 0^"^ in Eq. (0) to be a 
plane wave. Maxwell's equations (||) have the solution 

g-O'* = ,5^^;oe'('^=^-'=^) (23) 

for a slowly varying wave propagating in the x-direction with (we)^ — = 0. For a wave whose spatial variation 
is large compared to the size of the EIH surface surrounding the charge the surface integral equation (|^) yields the 
equation of motion 

mixUrer=-^qEoe'^'^ (24) 



which has the solution 

^l = l^ — ^Eoe^"'^. (25) 



It remains finally to compute the scattered field produced by the motion of our charge. One possible way to do 
this would be to introduce a model source term into the right hand side of the first of Eqs. (||) whose motion is 
characterized by x\ above. However, in keeping with the EIH philosophy of not specifically introducing sources into 
the field equations and because the EIH procedure makes it unnecessary, we will proceed by constructing a radiation 
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solution whose inner expansion matches on to the outer expansion of the near field (j) in Eq. ( |l7| ) with ta computed 
using x\ above. This outer expansion in inverse powers of r is given by 

= Yn 

(j,2 _|_ ^2 _ 2rri cos 6) 
= ^ + %riCos(? + 0(l/r3). (26) 



°+(^'(-^---» + ^(^/\r»^os. (27) 



The corresponding outer dipole field is given by 

a {W'\ 

er y er (er) 

whose inner expansion is 

a , Wjer) 

= ^ ( \i 

er yery 

where a and W[eT') are to be determined by matching. Comparing Eqs. (^9) and (|2l) we see that 



0(1) (28) 



a^eqx and W{eT) = eqin. (29) 

The gauge condition A'^,^= used to derive the wave equations for allows us to derive an expression for the 
vector potential part of corresponding to (j) given by 

A, = ^BlIl2fL, (30) 
r 

With the above expressions for and we can now calculate F^^, which in turn allows us to determine the 
Poynting vector 5''' = T"'' to be 

= {(1 + "3)m, (1 + m)n2,^{nl + nl)} (n,., f + 0(1). (31) 

where is a unit vector. The total integrated flux / is then calculated to be 

^oV- (32) 
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The Thompson scattering cross section a^, equal to / divided by the incoming flux Eq/Att, is thus given by 
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4^1 (33) 
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and is seen to differ from the usual expression for ctt by the factor of 1/R^. 

This result raises a number of complex questions. First and foremost, is it true or does some mechanism shield 
elementary processes such as Thompson scattering from the effects of the expanding universe? Do all cross sections, 
e.g. nuclear, experience this same effect? If the answer to both questions is yes, then a number of issues in stellar 
evolution and nucleosynthesis in the early universe will of necessity have to be reexamined. 
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